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We develop a framework for studying collective three-flavor neutrino oscillations based on the 
density matrix formalism. We show how techniques proven useful for collective two-flavor neu- 
trino oscillations such as corotating frames can be applied readily to three-flavor mixing. Applying 
two simple assumptions and the conservation of two "lepton numbers" we use this framework to 
demonstrate how the adiabatic/precession solution emerges. We illustrate with a numerical exam- 
ple how two stepwise spectral swaps appear naturally if the flavor evolution of the neutrino gas 
can be described by such a solution. For the special case where mu and tau flavor neutrinos are 
equally mixed and are produced with identical energy spectra and total numbers, we find that 
one of the spectral swaps in the three-flavor scenario agrees with that in the two-flavor scenario 
when appropriate mixing parameters are used. Using the corotating frame technique we show how 
the adiabatic/precession solution can obtain even in the presence of a dominant ordinary matter 
background. With this solution we can explain why neutrino spectral swapping can be sensitive to 
deviations from maximal 23-mixing when the "mu-tau" matter term is significant. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

It has long been recognized that, in addition to the 
conventional Mikheyev-Smirnov-Wolfcnstein (MSW) ef- 
fect [U, [3, Q , neutrino self-coupling can be important for 
neutrino flavor evolution when neutrino number densities 
are large I, i, i, 0, i, i, M M M MMM- Recently 
two-flavor neutrino oscillations in the core-collapse su- 
pernova environment have been intensively investigated 
[H [13, [S [H, US [m, m, m, m. These studies show 
that supernova neutrinos can indeed experience collective 
flavor evolution because of neutrino self-coupling, even 
when the neutrino self-coupling is subdominant com- 
pared to the MSW potential [iq . 

An important result is that collective two-flavor neu- 
trino oscillations can exhibit "stepwise spectral swaps" 
or "spectral splits" in the final neutrino energy spectra 
when the neutrino number density slowly decreases from 
a high value, where neutrinos experience s ync hronized 
oscillations [ij], towards zero (see, e.g., Ref. [lal)- When 
this occurs, u^'s appear to swap their energy spectra with 
Ux's at energies below or above (depending on the neu- 
trino mass hierarchy) a transition energy i5| x 2 (where the 
superscript "s" can stand either for "swapping point" or 
"splitting point"). Here Vx is some linear combination 
of V/^i and Vr- The phenomenon of spectral swapping is 
present in both the "single-angle approximation" , where 
flavor evolution along various neutrino trajectories is as- 
sumed to be the same as that along a representative tra- 
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jectory (e.g., the radial trajectory), and the "multi-angle 
approximation" , where flavor evolution along different 
trajectories is independently followed [13] ■ For the in- 
verted neutrino mass hierarchy, it also has been found 
that stepwise neutrino spectral swapping is essentially 
independent of the 2x2 effective vacuum mixing angle 
when this angle is small (see, e.g., Ref. [22l|). 

Collective two-flavor neutrino oscillations are best un- 
derstood with the help of neutrino flavor polarization 
vectors (3] or neutrino flavor isospins [16.]. Using the 
spin analogy, one can represent the flavor content of a 
neutrino mode by a spin vector in flavor isospace. In 
this analogy, the matter effect is described as a spin-field 
coupling, and the neutrino self-coupling is described as 
spin-spin coupling. The phenomenon of spectral swap- 
ping is a result of collective precession of all neutrino 
flavor isospins with a common angular velocity f2pr at 
any given neutrino number density [itI . [21[ . This collec- 
tive precession of neutrino flavor isospins is described the 
two- flavor adiabatic/precession solution [2S[2^. In this 
solution, all neutrino flavor isospins stay aligned or an- 
tialigned with a total effective field in a reference frame 
that rotates with angular velocity Qp-^. Numerical simu- 
lations have shown that, in the supernova environment, 
neutrinos can first experience collective MSW-like flavor 
transformation (in which the MSW effect is enhanced by 
neutrino self-coupling) and then subsequently the adia- 
batic/precession solution [13] • 

In the real world, however, there are three active neu- 
trino flavors. Some limited progress has been made 
on understanding collective three-flavor neutrino oscil- 
lations. The first fully-coupled simulation of three- 
flavor neutrino oscillations in the supernova environment 
showed a spectral swapping phenomenon similar to the 
two-flavor scenario except possibly with two swaps [28| . 
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Another simpHfied numerical study with a single neu- 
trino energy bin [29] reveals that collective neutrino os- 
cillations can be sensitive to deviations from maximal 
23-mixing when there is a dominant "mu-tau" term aris- 
ing from a high order contribution from virtual /i's and 
r's (e.g., Refs. 0, 0, EH). Ref. [32| has extended the 
neutrino flavor polarization vector notation to the three- 
flavor scenario and discussed the collective three-flavor 
oscillations as "factorization" of two two-flavor oscilla- 
tions. 

In this paper we develop a framework for studying 
collective three-flavor neutrino oscillations based on the 
density matrix formalism. Using this framework we flnd 
a generalized three-flavor adiabatic/precession solution 
and show how stepwise spectral swapping can appear as 
a natural result of such a solution. 

The rest of this paper is organized as follows. In Sec.HIl 
we develop a framework centered around a 3 x 3 reduced 
neutrino flavor density matrix. This density matrix is 
equivalent to an 8-component neutrino flavor vector, a 
generalization of neutrino flavor isospin and similar to 
the Bloch vector used in Ref. [l^]. We show how tech- 
niques important in studying collective two-flavor oscil- 
lations such as corotating frames can be applied in the 
framework. In Sec. IIIII we demonstrate how the three- 
flavor adiabatic/precession solution can be found using 
two simple assumptions and the conservation of two "lep- 
ton numbers" . We also illustrate with a numerical ex- 
ample how two spectral swaps can form from the adia- 
batic/precession solution when the total neutrino num- 
ber density vanishes. In Sec. lIVI we employ the corotating 
frame technique and show that the adiabatic/precession 
solution obtains even in the presence of a dominant mat- 
ter background. In particular, we show that, in the pres- 
ence of a large mu-tau term, neutrino spectral swapping 
becomes sensitive to deviations from maximal 23-mixing. 
In Sec. |V]we give our conclusions. 



II. EQUATIONS OF MOTION 

A. Density Matrix Description 

The flavor content of a neutrino (antineutrino) mode 
with momentum p is generally described by density ma- 
trix Pp {pp) 0. The diagonal elements of a density ma- 
trix are the occupation numbers of the neutrino eigen- 
states in a particular basis, and the off-diagonal elements 
contain the phase information relating to neutrino mix- 
ing. For a neutrino pure state described by flavor wave- 
function 




the density matrix is 

where is the neutrino number density in momen- 
tum mode p, and 01/1(1/2,1/3) are the amplitudes for the 
neutrino to be in the corresponding vacuum mass eigen- 
states. With this notation we have normalization 

E Kl' = i- (3) 

i=l,2,3 

The density matrix p for an antineutrino pure state is 
defined similarly. 

In this paper we will assume that the CP-violating 
phase is (5 = 0. A brief discussion of the effect of a non- 
vanishing CP phase is given in Sec. |Vl Because we are 
only interested in neutrino flavor transformation, we will 
assume that neutrinos are free streaming except for for- 
ward scattering on the background medium. For now we 
also will assume that there is no ordinary matter back- 
ground. For this case it is most convenient to work in the 
vacuum mass basis. This basis is implicitly adopted in all 
the following discussions except for Sec. lIVl where we will 
discuss the effects of the ordinary matter potential. To 
simplify the problem even further, we will assume that 
the neutrino gas is isotropic and uniform. This corre- 
sponds to the "single-angle approximation" in numerical 
simulations of ffavor oscillations of supernova neutrinos. 
It has been shown numerically [13, [IE [H, [ll] that both 
single-angle and multi-angle calculations produce similar 
neutrino spectral swaps. 

For an isotropic and uniform neutrino gas confined in 
a fixed volume, the equations of motion (e.o.m.) for neu- 
trino density matrix pp are 

iPp-[Hp,Pp]. (4) 
Here the Hamiltonian for neutrino mode p is 

Hp = H™^ + x/2Gf j (0(Pq - P;), (5) 

where Gp is the Fermi constant, and the vacuum term in 
the Hamiltonian is 

Hr = ^diag(m2,m2,m2) (6) 

with mf being the mass-squared eigenvalues correspond- 
ing to neutrino vacuum mass eigenstates For an- 
tineutrino density matrix pp one has 

iPp = [Hp,Pp], (7) 
with the Hamiltonian defined as 

Hp = H™^ + V2Gf j ^ (p-q - p* ). (8) 
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We note that Ref. has defined the antineutrino den- 
sity matrix as the complex conjugate of what one usually 
writes as a density matrix (i.e., p* — > pp). This notation 
leads to a slightly simpler version of the e.o.m. for both 
neutrinos and antineutrinos. However, for a vanishing 
CP phase, it is possible to treat neutrinos and antineu- 
trinos on an equal footing and combine Eqs. ^ and (171) 
into a single expression, as has been done in Ref. [l6| 
for the two-flavor mixing scenario. To see this, we note 
that a neutrino or antineutrino mode in an isotropic and 
uniform gas is completely characterized by 



= TTr(H™^A3) = ± 



At 



'21 



2E 



(9) 



where E ~ |p| is the energy of the neutrino or antineu- 
trino, upper (lower) signs are for neutrinos (antineutri- 



nos), Atoji 



mf is approximately the solar mass- 



squared difference Atoq , and A3 is one of the Gell-Mann 
matrices Aa (a = 1, ... ,8). Because the number density 
for a neutrino (antineutrino) mode z/p (Dp) is conserved 
for a neutrino gas in a fixed volume, we can define the 
total neutrino number density 



d^q 



Tr(pq + Pq), 



(27r)3 

and the normalized distribution function is 



dE 



Tr(pp) ifcc;>0, 
Tr(pp) if w < 0. 



(10) 



(11) 



Using Eqs. PH)) and pTjl we can express the integral in 
Eq. dH) as 



d^c 



:(Pq-Pq) 



(12) 



In Eq. p2)) we have defined the "reduced neutrino flavor 
density matrix" for neutrino mode u: 



Pq if > 0, 
-p* if < 



which has normalization 



Tr(£i^) 



-1 if w > 0, 
-1 if < 0. 



(13) 



(14) 



Using Eqs. (j4|), ([7| and ([T2]) we find the e.o.m. for Qu;-. 

'^Q<^ = [Hc^, Quj]- (15) 
The Hamiltonian for neutrino mode u is 

- H^^ + ^igtot, (16) 
where in the vacuum mass basis 



2 V3 



(17) 



Here we define 



V2G 



(18) 



This parameter dictates the strength of neutrino self- 
coupling. The total neutrino flavor density matrix is de- 
fined to be 



Qtot 



(19) 



In Eq. pB]) we have left out the trace term for 
(which is irrelevant for neutrino oscillations), and we have 
defined the oscillation parameter k to be 



zp^Tr(H™'= 



= ±- 



1 

'2E 



mo 



As), 
{ml 



(20a) 
(20b) 



where the upper (lower) signs are for neutrinos (antineu- 
trinos). Because the atmospheric mass-squared differ- 



i^tm is much larger than Attjq, one has 



(m 



(21) 



where the plus (minus) sign is for the normal (inverted) 
neutrino mass hierarchy. 

We note that does not change with time if there is 
no inelastic scattering of neutrinos. We also note that 
Eq. (fT^l) is actually more generally valid than Eqs. @ 
and ^ so long as the neutrino gas stays isotropic and 
uniform. For example, this would be true for an homo- 
geneous and isotropic early universe, i.e., the Friedman 
solution. The "single-angle approximation" (see, e.g., 
Ref. [17j) for fiavor evolution of supernova neutrinos is 
essentially equivalent to this scenario. For this case, the 
fiavor content of a neutrino propagating along any trajec- 
tory at a given radius is assumed to be identical to that 
of a neutrino with the same energy propagating along a 
radial trajectory at the same radius. In this approxima- 
tion the flavor evolution of neutrinos as a function of time 
t can be represented as the flavor evolution of neutrinos 
propagating along the radial trajectory as a function of 
radius r. In addition, one can define the effective total 
neutrino number density at each radius as 



D{r/R,) ^ 



dEUiE). 



(22) 



This takes account of both the geometric dilution and 
(partly) the anisotropy of supernova neutrinos. In 
Eq. ((22|l Ri, is the radius of the neutrino sphere, the 
geometric factor is 



D{0 = l{l-Vl^)^ 



(23) 



Li,, (E^) and /i/(£') are the neutrino luminosity, average 
energy and normalized energy distribution function for 
species v at the neutrino sphere, respectively, and the 
summation is over all neutrino species (including both 
neutrinos and antineutrinos). 
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B. Flavor Vector Description 

The flavor polarization vector P Q and neutrino flavor 
isospin [16| are important techniques that have been 
used extensively to describe the flavor mixing of neutri- 
nos in two-flavor mixing scenarios. These notations have 
helped in visualizing and giving insights into the problem 
of collective neutrino oscillations. 

To generalize these notations to the 3x3 case, we note 
that a 3 X 3 Hermitian matrix A can be expressed as the 
linear combination of the identity matrix I and Gell-Mann 
matrices Xa- 



where 



A=iTr(A)l+^A('^)^, 



EE Tr(AAa) 



(24) 



(25) 



can be viewed as the a'th component of an 8-dimensional 
vector A. In particular, "flavor vector" 



g _ (Jl) „(8)nT 



(26) 



is the generalized version of neutrino flavor isospin s^^.^ 

Because Tr(pt»j) is flxed by the normalization condition 
in Eq. p4|) . flavor vector is fully equivalent to the 
density matrix g^^. In particular, the number densities of 
vi , V2 and 1^3 in mode uj can be expressed in terms of q\ 
and Q, 



(3) 



(8). 



(27a) 
(27b) 



n.3(c.)=<°V.Q-^^>(f^). (27c) 
Noting the difference between the definition of for 



^ Ref. I32I appears while this manuscript was in preparation. 
Ref. I32II has proposed a three-flavor version of the Bloch vec- 
tor which is a generalization of the two-flavor polarization vector 
defined in Ref. Q . The difference between the three-flavor Bloch 
vector and the flavor vector defined here is similar to that be- 
tween the two-fiavor polarization vector and the neutrino flavor 
isospin defined in Ref. ITsll . In the fiavor vector description, the 
directions of flavor vectors and flavor isospins for antineutrinos 
are intentionally reversed so that the e.o.m. for flavor vectors 
and flavor isospins (for both neutrinos and antineutrinos) can be 
written within a single expression. Likewise, the density matrix 
defined in Ref. |32( is different from the flavor density matrix 
in this paper by a sign for antineutrinos. This is of course a 
notation difference and does not affect the physical results. 



neutrinos and antineutrinos [Eq. ((T3)) ]. we have 

1 



1 



2V3 



,(8) 



(28a) 
(28b) 
(28c) 



The relation between uj and the energy of a neutrino or 
antineutrino is described in Eq. ([9]). 

One can define the cross and dot products of two 8- 
dimensional vectors A and B to be 

A X B = -i Tr([A, B] Aa)e, = fabcA'^^^B'^^^Ba, (29a) 

A • B = 2Tr(AB) - ^Tr(A)Tr(B) = A^^^B^^^ (29b) 

where is the a'th unit vector in the 8-dimensional fia- 
vor space, and fabc are the antisymmetric structure con- 
stants of SU(3): 



Aa Af, 

Y'Y 



i/ah 



A. 



(30) 



The summation over Gell-Mann indices are implicitly as- 
sumed in the above equations when they appear twice in 
the subscripts or subscripts. Using Eqs. (|15p and (|29ap 
we can write the e.o.m. for fiavor vector as 



at 



2k , 



(31a) 

/OO 
-00 

(31b) 



The first term in Eq. (j31bp corresponds the precession 
of around an "external field" , and the second term 
corresponds a "spin-spin anti-coupling" with strength /z. 

Eq. (|31ap makes it clear that the change in g,^ is or- 
thogonal to itself, and, therefore, the magnitude of 
this quantity does not changed. We have 



le^P = Xl[Tr(ec^Aa)]^ = const. 



(32) 



This is a natural result that also follows from the fact 
that Eqs. ^ and (JTlmaintain the coherence of pp and 
Pp. Following Ref. [l6j we define the effective energy of 
the system to be 



- I duj duj' fuifui'Qu, ■ Qu. 



(33) 



Clearly the effective energy £ of the system is conserved 
if 71*°* is constant. 
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Although flavor vector seems to behave in a way 
similar to flavor isospin s^^ in two-flavor mixing scenar- 
ios, there are fundamental difl^erences between the 8- 
dimensional flavor vector space and the 3-dimensional 
flavor isospace. For example, two 8-dimensional vectors 
A and B are "perpendicular" or "parallel" to each other 
if A-B = or A X B = 0, respectively. Because there are 
two linearly independent generators of SU(3) that com- 
mute with each other, for any vector A one can always 
find another vector A' which is both "perpendicular" and 
"parallel" to A. Consequently, generally speaking 



B ^ A 



A B 



(34) 



even if B is "parallel" to A. 

The existence of two linearly independent and com- 
muting generators of SU(3) has another important con- 
sequence. Because [A3,A8] = 0, rotations around 63 and 
eg can be viewed as independent. In particular, the first 
term in Eq. (j31bp can be interpreted as simultaneous but 
independent precession of flavor vector q,^ around 63 and 
eg with generally different angular velocities. 

Although the density matrix and flavor vector descrip- 
tions are equivalent, the rotation of a flavor vector in 
the 8-dimensional flavor space is not easily visualizable. 
Therefore, we will base our discussions mostly on the 
density matrix formalism with intermittent references to 
the flavor vector notation where it seems convenient. 



C. Conserved "Lepton Numbers" 

Multiplying Eq. (jlSp by f^j and integrating it over lj 
we obtain the e.o.m. for (Stot^ 

mot = 7; dw/c^a)[£»i^, A3] + — = / duj f^K[g^,Xs]. 
^ J -00 V o J -00 

(35) 

Eq. (j35p is not a closed equation from which we could 
solve for Qtot- However, because A3 and Ag commute with 
each other, it is clear that the two "lepton numbers" 



■^3 = Qi 



(3) 



and 



(8) 



dw /^Tr(£)^A3) (36a) 



dw/„Tr(gi„Ag) (36b) 



are constants of the motion. Because Tr(£itot) does not 
change with time, the lepton number (fraction) in each 
vacuum mass eigenstate is individually conserved. 



D. Corotating Frame 

In the density matrix description, changing from the 
static frame to a corotating frame corresponds to a trans- 



formation 

Quj ^ Quj = e p^e , [61 ) 

where Hcor is a 3 x 3 Hermitian matrix that is common to 
all neutrino modes and does not change with time. The 
flavor density matrix in the corotating frame satisfies 
an e.o.m. similar to that in Eq. p^ : 



ipw — [Hcj, Qij], 



where 



-,-iHc 



(38) 



(39) 



is the Hamiltonian for neutrino mode uj in the corotating 
frame associated with Hcor. 

A very special set of corotating frames corresponds to 
simultaneous rotations around §3 and eg with angular 
velocities 51 and 2K/\/3, respectively, and 



n— - K-^. 
2 Vs 



(40) 



Because H™'^, A3 and Ag commute with each other, H^^ 
in these special corotating frames takes the same form as 
in Eq. except for the replacements 

u} ^ Lu - ft, (41a) 

K-^ K~K, (41b) 

Qtot ^ Ptot = e'^-*gtote-'"-'*. (41c) 

Also, because the occupation numbers of each vacuum 

eigenstate in a neutrino mode lu are determined only by 

(3) (s) 

Qui and Quj , the probability for the neutrino mode to 
be in the i'th vacuum mass eigenstate in these special 
corotating frames, |a^;(a')p, is the same as the probabil- 
ity for the neutrino mode to be in the same eigenstate 
in the static frame, |a^;(a')p. Therefore, lepton num- 
bers L3 and Lg are not changed by the corotating-frame 
transformation either. 



III. ADIABATIC/PRECESSION SOLUTION 
AND STEPWISE SPECTRAL SWAPPING 

We now seek a natural extension to the adia- 
batic/precession solution presented in Ref. [2^. This so- 
lution, like the 2x2 case, will be essentially a quasi-static 
solution that, for given n*°', is the same as the "static" 
solution which satisfies the precession ansatz. There is 
a family of static solutions for each value of n*°*. So 
long as changes slowly, a particular solution in each 
family parametrized by n^°^ is uniquely determined by 
the initial conditions and by the adiabatic ansatz. This 
gives the adiabatic/precession solution. We will discuss 
the precession ansatz and three-flavor synchronization in 
Sec. nil Al In Sec. IIIIBI we will discuss the adiabatic 
ansatz and outline a formal procedure for obtaining a 
three-fiavor adiabatic/precession solution. In Sec. IIII CI 
we will illustrate with a numerical example how step- 
wise spectral swapping arises from the three-flavor adia- 
batic/precession solution as n*°' 0. 
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A. The Precession Ansatz and Synchronization 

The precession ansatz is that, for constant it is 
possible to find a Hermitian matrix 



V3 



(42) 



such that the flavor density matrix is static in the 
corotating frame associated with Hcor, i-e. 



(43) 



From the arguments in Sec. Ill Dl if the precession ansatz 
is satisfied, ja^. (a;)p does not change with time. In this 
sense, a solution that satisfies the precession ansatz is 
"static" . 

When the adiabatic ansatz is satisfied, each flavor vec- 
tor will precess uniformly around 63 and eg. In other 
words, the system is in a state that is symmetric about 
63 and eg. Because the e.o.m. for [Eqs. (PT|) ] pos- 
sess the same symmetry around the 63 and eg axes, we 
expect "static" precession solutions to exist for constant 
n*°* with appropriate initial conditions. In particular, 
this symmetry obtains approximately in dense neutrino 
gases where 



%/2Gf<°* > |k| > \uj\ 



(44) 



for most neutrino modes. This conclusion can be shown 
as follows. 

When Eq. is satisfied, the conserved effective en- 
ergy of the system is [Eq. 



and, therefore, 



I Qtot I 



const. 



(45) 



(46) 

Ptot and Lg = g'^^ 



Because the two lepton numbers L3 
are conserved (Sec. Ill Ct . the total fiavor vector ^tot must 
precess simultaneously around 63 and eg according to 



Qtot X ( 17^63 + — ^1 



Also in this limit, Eq. ((3T|) becomes 
d 



-P-Quj X ^tot- 



(47) 



(48) 



Eqs. (|47|) and (|48)) suggests a simple geometric picture 
for the flavor evolution in dense neutrino gases. On short 
time scales (At ~ M^^), each flavor vector g^^ precess 
rapidly around the total fiavor vector ^tot- On large time 
scales (At ~ jfi 



pr I 



), all flavor vectors precess 
slowly around 63 and eg with angular velocities and 

2 AT^ / \/3, respectively. This is analogous to the synchro- 
nization phenomenon in two-fiavor mixing scenarios . 



B. The Adiabatic Ansatz and 
Adiabatic/Precession Solutions 



When the precession ansatz in Eq. (|43|) is satisfied, it is 
possible to find a unitary matrix X^^ that simultaneously 
diagonalizes both and g^: 

XcjHi^X^^ = diag((I'L,<I'M,'^H), 



Xcj^c^X^^ = ±diag(|c 



(49a) 

laPMlMa^^Hp), (49b) 

where wl < i^M < i^H are the eigenvalues corresponding 
to the eigenstates of H^^, \i'^{ui)) and |i/H(u;)). 

In Eq. (j49b|) the plus (minus) sign is for neutrinos (an- 
tineutrinos) . The adiabatic ansatz is simply that 

jaj;^ (a;)|2 ^ const. (/ = L, M, H) (50) 

as slowly varies with time. 

As discussed in Ref. [l^], this "adiabaticity" criterion 
connects neutrino systems in different corotating frames 
at different values of n*°* . Note that this adiabaticity cri- 
terion is different from the meaning of adiabaticity usu- 
ally adopted in the literature, e.g., when discussing the 
MSW mechanism, which is always based on the static 
frame. Following Ref. [2^, we argue here that the adia- 
batic ansatz can be satisfied if, for each neutrino mode 

rotates at a speed much slower than precession rate 
of gi^ around H^^, i.e. 



7 ^ 



|H^ X dH^/dt| 



< 1. 



(51) 



The adiabatic/precession solution can be obtained for- 
mally by employing the following procedure: 

1. At any n*°* find for each neutrino mode uj a uni- 
tary matrix X„ that diagonalizes H^^. Matrix X^^ is 
expressed as a function of (gtot, f^pr, ATpr). 

2. For given initial values of |a5j(u;)p, find g^ = 



±Xidiag(|< 



'£>LlMasMlM«i>Hl )Xw as a function of 



(^tot, ^^pr, A'pi-). 

3. For given initial lepton numbers L3 and Lg find 
g^tJt — ^3 and g[fl = Lg. From the definition ptot = 
J^^dLuf^g^{gtot,^pr, ATpr), solve for the precession 
angular velocities f2pr and ATpr and the remaining 
components of ^tot- 

4. Find g^j for each neutrino mode uj using the expres- 
sion gui{gtot,^pn Kpr) obtained in stepO 

We note that the above procedure actually gives a set 
of equivalent solutions. This is because the precession so- 
lution is symmetric around 63 and eg. If (ptot, ^^pr, A'pi) 
is a solution, {g'toti ^pn ^pr) is also a solution, where g't^t 
is related to ptot by two arbitrary phases (f)^ and (/)g: 

. A, 



^tot 



exp [ -103 Y ^ Y ) ^tot exp M'/'s y + i^s y 



One can fix these two phases by, e.g., choosing gtol 

~{7) 

Qtot 



(52) 



0. 
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C. Stepwise Spectral Swapping 



where 



Neutrino flavor mixing becomes very simple in the adi- 
abatic/precession solution presented in Sec. IIII Bl when 
n*°* 0. In this limit is diagonal in the vacuum 
mass basis: 



Ha;|„tot^o = diag((2;i,a)2,cj3)„tot^o, 



where 



-0 = -{n-K^,), 



(53) 

(54a) 
(54b) 
(54c) 



with Slpj and K^^. being the collective precession angular 

velocities as n*°* 0. Eq. ()53p shows that has 
3 critical values of lo at which any two of its eigenvalues 
are equal. These critical values are 



""'31 



'32 



K° 

pr 



pr 



2"pr 
2"pr 



(55a) 
(55b) 

(55c) 



In practice, however, two of the critical points are usually 
indistinguishable 



^0 = "pr; 



Wo 



OJn 



±- 



Ar 



I ^'^pri 



(56a) 
(56b) 



because Am^^^ ^ Attzq and, therefore, \Kpj.\ \^pr\- 
In Eq. (|56b|) the positive (negative) sign is for the nor- 
mal (inverted) neutrino mass hierarchy. A critical value 
Lu^ corresponds to the energy where "stepwise spectral 
swapping" or a "spectral split" occurs. If a;** > 0, the 
spectral swap occurs at neutrino energy 



At 



(57) 



in the neutrino sector. If cj^ < 0, the swap is located in 
the antineutrino sector at energy E^. 

We can illustrate stepwise spectral swapping by us- 
ing the following test example. We assume a bare, hot, 
spherical neutron-star that isotropically emits neutrinos 
directly into the vacuum from its infinitely thin neutrino 
sphere. We adopt the single-angle approximation and 
take the radius of the neutrino sphere to be Ri, = 30 
km, and the luminosity for each neutrino species to be 
= 10^^ erg/s. We take the energy spectra of neutrinos 
to be of the Fermi-Dirac form: 



ME) 



1 



T3F2(7y,)exp(S/r„-77,) + r 



(58) 



Fkiv) 



X dx 



exp(a; — t]) + 1 



(59) 



We take degeneracy parameters to be the same for all 
neutrino species, ry^ = 3, and we choose to be such 
that the average energies for various neutrino species are 
{E,J = 11 MeV, {E^J = 16 MeV and {E,J = (iJpJ = 
{Ei,^) = (Efj^) = 25 MeV, respectively. For the neutrino 
mixing parameters (see, e.g., Ref. [33| for our conven- 
tions) we take 612 = 0.6, 6*13 = 0.1, 6*23 = 7r/4, ^ = 0, 



-3 X 10-3 eV^ (in- 



Am|i = 8 X 10-5 eV^ and Aml^ 
verted neutrino mass hierarchy). 

We define P(i^a Vi) and P(j?Q Di) as the proba- 
bilities for neutrinos and antineutrinos that are initially 
in pure a flavor state at the neutrino sphere to end up in 
the i'th vacuum mass eigenstate as n^°^ — > 0. In Fig. [1] 
we show neutrino conversion probabilities P{i'a ~* ^i) 
and P{va — > i^i) as functions of neutrino energies in our 
test example. We observe that while P{va Vi) shows 
two swaps at E^y ~ 5.2 MeV and El^^ ~ 8.4 MeV, re- 
spectively, P{Da Vi) shows no swap at all. This phe- 
nomenon can be explained using the adiabatic/precession 
solution discussed above. 



TABLE I: The correspondence between neutrino states \i'i{ijj)) 
and \va) or \h'oi) at r — R^, for the inverted neutrino mass 
hierarchy case with mixing angles 613 ~ and 623 ~ 7r/4. 





uj <0 


0<u;< culZ^ 


^ sync 


\M^)) 

\i>M{tLi)) 
l'^(^)) 


^i\v,) - \v.)) 

^i\v^) + \v-r)) 
\Vc) 


i=,i\v,) + \v.}) 

\Vc) 


i=,iK) + \v.}) 

\Vc) 



TABLE II: The correspondence between neutrino states 
\i}i(Lo)) and \vi) or Ivi) as nl;°* — > for the inverted mass hier- 
archy case and for swap points with hierarchy luq > tjlt^ > 0. 





LO<0 


<UJ < Ultn, 






\M^)) 


\V2) 


\V2) 




\y3) 


\i>M{iL!)) 


\vi) 


W) 


\V2) 




|wi(^)) 


\V3) 


l^3> 




\V2) 



TABLE III: Nonzero neutrino conversion probabilities 
P(ya Vi) and Piva Vi) based on Tables HI and HTl and the 
adiabaticity ansatz in Eq. (|66|l . 



energy range 


nonzero conversion 


probabilities 


any 


P{ve -> V3) = 2P{u^ 


T —> Vl,2) = 1 


El, > -EJtm 


P{Ve V3) = 2P{l^^ 


T Vi,2) = 1 


Eq < Ev < -EJtm 


P{Ve Vx) = 2P(l/^ 


T V2:3) = 1 


E,j < Eq 


P{Ve V2) = 2P{l^^ 


T Vi,'i) = 1 




FIG. 1: (Color online) Conversion probabilities P{i'c, — » Vi) (left panels) and P{i'a i'i) (right panels) as functions of neutrino 
and antineutrino energies and respectively, in the bare, hot neutron star example. The top, middle and bottom panels 
are for neutrinos that are initially in pure e, /i and r flavors, respectively. The solid, dashed and dot-dashed lines are for 
neutrinos that end up in the 1st, 2nd and 3rd vacuum mass eigenstates, respectively, when n^°^ 0. 



Assuming that the precession ansatz in Eq. is sat- 
isfied, we can diagonalize the Hamiltonian H(j|r=fl„ in the 
corotating frame in which aU flavor vectors are static. We 
denote wl < oJu < as eigenvalues corresponding to 
the eigenstates ji'L), \vm) and li'H) of H^|r=_R„, respec- 
tively. Because is very large at the neutrino sphere, 
we have 

\/2GfK, -npj > |Lj-n~|,|At-if~| (60) 
for most neutrino modes. Therefore, 



|i>H(w))r=i 



l^e) if W > 0, 

iPp) if w < 0. 



(61) 



We note that fip^ and Kp,- are essentially a kind of aver- 
age of Lo and k in the system. In our example, neutrinos 
(instead of antineutrinos) are the dominant species and 
the neutrino mass hierarchy is inverted. So we expect 



fipr > > at any value of We can diagonalize 
the ^r-submatrix of H™'^ in the flavor basis and obtain 



\i>L{^))r=R^ - -^(K) - Wr)), (62a) 



(62b) 



a u! < 0, and 



1 

1 



\l'L{(^))r=B„ ^ —iWt^) T Wt)), (63a) 



\i>Mil^))r=R^ ^ —(Wfi) ± I^^t)), (63b) 



if a; > 0. In Eq. (j63p . the upper and lower signs are for 
the cases where uj is smaller or larger than 



, sync ^ T^oo 

^atm ^^pr \ 2 ' 



(64) 
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respectively. Here |£'L(w))r=i?,„ and \vMi^))r=R^ are 
"equal mixes" of ji'^) and Ij^t-) or 1/?^) and These 
results are summarized in Table HI 

Far away from the neutron star where n^°^ — > 0, we 
obtain |i>;)„i.ot_^o using Eqs. ([53)1 and ((54)) . For a; < 
(antineutrinos) we have 

|z>L(w))<ot^o = \i^2), (65a) 
\i>u{uj)) ni^'-^Q = (65b) 
|^'H(w))„t_ot^O = |i^3>- (65c) 

Assuming that uJq > i^l^^-^ > 0, we summarize the corre- 
spondence between p/(cj)) and for n|;°* — > in Table 

m 

If the adiabatic ansatz in Eq. (f50|) is satisfied, we have 

|a£>,(t^)lr=_R„ = |a£>,(^)lr»t°'^o (66) 

I 



for each neutrino mode u. Using Eq. (|66p and Tables 
U and [ni we can obtain the neutrino conversion prob- 
abilities Pii'a i^i) and P{Pa — i^i) for each neutrino 
mode uj. In Table [TTTl we summarize the values of nonzero 
neutrino conversion probability when the flavor evolution 
of the neutrino gas follows the adiabatic/precession solu- 
tion. The results in Table [TTTl and Fig.[l]are in good agree- 
ment. The neutrino spectral swapping feature in Fig. [1] 
indeed can be explained by the adiabatic/precession so- 
lution. 



As in two-flavor scenarios, the swapping energies Eq 
and i?atm <^an be obtained by using the conservation of 
lepton numbers and by assuming that the spectral swaps 
are infinitely sharp in energy. From conservation of = 

Ptot we have 



r=R„ 



dE {2U^ - 2UJ - / dE {2U^ - 2U, 



(67a) 
(67b) 



where fv{E) is a distribution function which satisfies 

f,[E) cx f,{E) (68) 
with normalization condition 



V rdEU{E)^i, 

„ Jo 



(69) 



and 



L = U = h^u. = u.- (70) 

From Eq. ([57)1 we can determine that i?atm — ^-4 MeV. 
We also can obtain Eq from the conservation of L3 = 



(3). 

Stof 



d£;cos20i2(/,, -/pj, (71a) 



/-^s 



dE{U^-U^ 



El, 



dE{U^-UJ. 



(71b) 



This implies E^q ~ 5.3 MeV. The values for El^^ and 
Eq derived from the adiabatic/precession solution are 
also in good agreement with the numerical results shown 
in Fig. dl 

We note that, in this numerical example, neutrinos of 
the fi and r flavors are equally mixed and have identi- 
cal energy spectra initially. In this particular case, the 



two-flavor approximation with effective mixing parame- 



ters Ar 



-^™atm and 



013 produce a similar 



spectral swapping feature around i?|x2 — -^atm- In the 
two- flavor mixing scenario, the value for -E|x2 also de- 
termined by the conservation of a lepton number 



^2x2 — J^l + h, 



I1-I2, 



where 



dEiu^^UJ, 



(72a) 
(72b) 

(73a) 



l2= dE i,U^ - ,UJ - / dE {U^ - /,J. (73b) 

From Eq. ([7^ it is easy to see that the integral I2 = 0- 
Comparing Eqs. ([67)1 and ([72)) . we see that a solution 
for -^2x2 ill Eq. ([7^ is a solution for i?atm — -^2x2 
Eq. ([57)) . Therefore, the two-flavor spectral swap phe- 
nomenon is completely consistent with the three-flavor 
spectral swap phenomenon at the atmospheric mass- 
squared-difference scale. 

We also note that, in this example, there would ex- 
ist only one spectral swap at energy E^, ~ El^^ if 
we had chosen the other hierarchy for swap points, i.e. 
'^atm > '^Q > 0. It is generally not possible for a single 
spectral swap to satisfy the conservation of both lepton 
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numbers and, therefore, this hierarchy for swap points is 
not physical. 

The example given here can be generalized into a 
generic procedure to predict neutrino spectral swaps: 

1. Diagonalize H^^ when is large, and find the cor- 
respondence between \i'i{uj)) and \va) or as was 
done in Table U 

2. Estimate the approximate locations of the swap 
points as nj^°* — > 0. The swap points are expected 
to be in the neutrino sector if the system is domi- 
nated by neutrinos instead of antineutrinos. Pick a 
hierarchy for swap points, i.e., whether E^^^ > Eq 

or i^Itm < ^0- 

3. Obtain from Eq. (|54p the correspondence between 

and {vi) or for n*°* — > 0, as was done in 
Table m 

4. Use the results in step [T] and [3] and the adiabatic 
ansatz in Eq. (|66p to find neutrino conversion prob- 
abilities P{h'a Vi) and P{va ^ i^i), as was done 
in Table Uni 

5. Find lepton numbers L3 and Lg as n'°* 0. These 
will be functions of Eq and -E^tm when the initial 
energy spectra fi, (E) and the results in step S] are 
used. 



becomes 



6. Solve for E% and E: 





conservation L 



3\t= 



^atm by using lepton number 
=0 = -^3Utot^o(£'0,£^atm) and 
i8|t=o = ^8|<°t^o(£'0j -E^Itm)- If no consistent so- 
lution can be found, pick the other hierarchy for 
swap points -Eq and E|jj,j in step [2] and repeat steps 



(74) 



where H™''** is the Hamiltonian contribution arising from 
neutrino-electron forward scattering. In general this 
term will give different refractive indices for Vp jvf^ and 
v^^t/v^^t P, Ignoring the trace term, in the fiavor 

basis we can write 



,_|matt ^ ^G'p„^diag(re,0,0), 



(75) 



where rib is the number density of baryons, and Yg is 
the electron fraction. MSW resonances can occur if rig — 
UhYe is small and comparable to \uj\/^/2Gf or \k\/V2Gf. 
In this section, however, we will assume rig to be constant 
and very large for most neutrino modes: 



(76) 



Because H™*^** is independent of uj, it vanishes in the 
corotating frame picked out by 



IV. ADIABATIC/PRECESSION SOLUTIONS IN 
A DOMINANT MATTER BACKGROUND 

A. Effects of Neutrino-Electron Forward Scattering 



(77) 



In the presence of ordinary matter, Eq. (|15p is still 
valid except that the Hamiltonian for neutrino mode lu 



In this corotating frame, the vacuum Hamiltonian for 
neutrino mode uj becomes 



/r 



2„2 



At 



ml 



hi2{t) 

2 2 



hiz{t) 



+ 777ICI2) \{ 



(78) 



\{Tnl + mlsl2 + 7nlcl2)) 



in the fiavor basis, where Cy = cosSy, ~ sin^y , and ^13 ~ and 623, — 7r/4. Because 

hi2(t) and hisit) are functions that oscillate with angular _ _ _ 

frequency V^Gpne. In deriving Eq. ^ we have taken V2GFne > I^i2(i)|, \hi3{t)\, (79) 

we expect hi2{t) and hi^{t) to average to and, therefore, 
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to have little effect on neutrino flavor evolution. 
Setting hi2{t) and hi^{t) to 0, we can diagonalize 



UJ 



diag(TOf 



where 



2 2 
2 



2 2 
■^12^21 



2,2 2 
12'™! + C]^2'^2i 



(80) 

(81a) 
(81b) 
(81c) 



are the effective mass-squared values for neutrino states 



1 



I^r)), 
Wr)). 



(82a) 
(82b) 

(82c) 



We note that the flavor density matrix in the coro- 
tating frame associated with Hcor — H™''** obeys an 
e.o.m. similar to that obeyed by g^^ in vacuum [Eq. (fT5|) ]. 
except for small perturbations occurring on very short 



time scales {At '--^ Gp 



The only difference is that 



the presence of a large net electron background breaks 
the "degeneracy" between Ve/i^e and v^^^T-jv^j^^T. There- 
fore, we can obtain adiabatic/precession solutions using 
the same procedure listed in Sec. IIIIBI but with the re- 
placements 



and TO? 



(83) 



Likewise, the conserved lepton numbers should be calcu- 
lated in the basis instead of the \vi) basis. 



B. Effects of Virtual Charged Leptons 

At very large matter density, virtual \i and r states 
contribute to a higher order correction to neutrino re- 
fractive indices, and H'°^" is found to be H H [Ml 



nbdiag(ye, 0, Fr) 



(84) 



in the flavor basis, where n\Xr gives the effective net r 
lepton abundance. If the matter density is so large that 



(85) 



we can again employ the corotating frame as in Eq. (j77p . 
Ignoring the rapidly-oscillating off-diagonal elements di- 
agonalizes Hc^ in the flavor basis in this corotating frame: 



where 



Ato^i 



diag(TOf ,TOf,TOf), 



112 2 / 2 2,2 2\|2 2 
"^2 -C23(Sl2mi+Ci2TO2) + S23'^3> 



112 2 / 2 2,2 2\ 



2 2 
^23 "^3 



(86) 

(87a) 
(87b) 
(87c) 



are the effective mass-squared value for the neutrino 

states 

10 ^ We). (88a) 

Wi) ^ K), (88b) 



Therefore, the adiabatic/precession solution still obtains 
in the presence of both large electron and large effective 
tau abundances. In this case, however, the degeneracy 
among Vf^jv^^ ^til^ti and Vt/i^t is completely broken. In 
this case, the adiabatic/precession solution is best ob- 
tained in the flavor basis. The conserved lepton numbers 
should also be calculated in the flavor basis. 

Because the absolute neutrino masses are irrelevant for 
neutrino oscillations, in the inverted neutrino mass hier- 
archy case we can take mi ~ to| ~ Ato^^j^ > TO3 = 0. 
When nhYr is large, from Eq. ([57)1 it can be seen that the 
mass-squared eigenvalue for ji'") is always heavier than 
the mass-squared eigenvalues for m') and ji/g ), while the 
effective 23-mass-hierarchy, or the sign of 



Am: 



112 _ 

'32 



112 

Too — Too 



//2 



(89) 



depends on whether ^23 is larger or smaller than 7r/4. 
Similarly, for the normal neutrino mass hierarchy case, 
we can take to§ ~ Ato^^j^^ > TOi ~ TO2 — 0, and again 
the effective 23-mass-hierarchy depends on whether ^23 is 
larger or smaller than 7r/4, although in a reversed fashion. 
According to the discussion in Sec. IIII C[ the flnal energy 
spectra for v'^jv'^ and v'^^jv'l (and, therefore, the spectra 
for Vfi/vi^ and v-t/v-t) interchange with each other when 
6*23 rotates from the first octant to the second octant. 
In other words, the final neutrino energy spectra can be 
sensitive to deviations from maximal 23-mixing. This can 
be illustrated using the toy model discussed in Ref. [29l |. 

In this toy model the neutron star emits only v^. and 
i/g with the same energy into a thick matter envelope 
where rtbFr is large. In this model it is assumed also 
that nii^jriy^ = 1 -I- e at the neutrino sphere with e > 0. 
If the neutrino gas follows the adiabatic/precession solu- 
tion, then P(z^e v'l) and PiVf^ v'l) must be either 
or 1 except at swapping points (see Sec. IIII C|) . However, 
it is clear from the conservation of lepton numbers in the 
flavor basis that v^, 's can not have been fully converted to 
other flavors. Because neutrino (instead of antineutrino) 
is the dominant species in the system, the spectral swaps 
occur in the neutrino sector. Therefore, Pe's are com- 
pletely converted into antineutrinos of another flavor as 
n*°* — > 0. In this case, the precession ansatz in Eq. (|43|) 
is trivially satisfled as n*°* ^ in the corotating frame 
defined by 



Hvac 



(90) 



where W^'^ is the vacuum term in the Hamiltonian for 
the neutrino mode with energy Ey. In other words, both 
swapping points have collapsed into one which is located 
in the neutrino sector and at energy E^ = Ey. Following 
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the discussions in Sec. IIII CI we can use the conservation 
of lepton numbers in the flavor basis to find that the 
neutrino conversion probabilities are 



t^yve - 




e 

1 + e' 








-0, 


(91b) 






1 

1 + e' 


(91c) 






= 0, 


(91d) 


P{i>e- 




-0, 


(91e) 


P{l>e- 




= 1 


(91f) 



for the inverted neutrino mass hierarchy and ^23 < 7r/4. 
For the inverted neutrino mass hierarchy and 023 > 7r/4 
we can obtain neutrino conversion probabilities which 
are similar to those in Eq. (|9ip but with <-!■ and 
i^fi Vr- This is exactly what has been observed in the 
numerical calculations for the toy model at r ~ 400 km 
(Fig. 2 in Ref. ^29,]). 

V. CONCLUSIONS 

We have developed a framework for studying collective 
three-flavor neutrino oscillations. Important techniques 
in studying collective two-flavor oscillations such as coro- 
tating frames can be applied readily to three-flavor sce- 
narios in this framework. We have shown that the three- 
flavor adiabatic/precession solution obtains when both 
the precession ansatz and the adiabatic ansatz are satis- 
fied. If the flavor evolution of a neutrino gas is described 
by the adiabatic/precession solution, the final neutrino 
energy spectra will exhibit the stepwise swapping phe- 
nomenon. We have shown that stepwise spectral swap- 
ping appears in a numerical example in which neutri- 
nos are directly emitted from the neutrino sphere of a 
bare neutron star into vacuum. For this special exam- 
ple, because the neutrinos in /i and t flavors are equally 
mixed and have identical energy spectra initially, the adi- 
abatic/precession solutions for both the three-flavor and 
the two-flavor scenarios produce the same spectral swap- 
ping at the atmospheric neutrino mass-squared-difference 
scale. In more general cases, however, this may not be 
the case, and the full 3x3 mixing framework should be 
employed. 

Strictly speaking, the adiabatic/precession solution ob- 
tains only when the CP-violating phase 5 = 0. This is 



because, if (5 ^ 0, the unitary transformation matrix U 
connecting the flavor states and vacuum mass eigenstates 
is not real, and 

UplU^ ^ {UppU^r and Up^U^ ^ {UppU^)*. (92) 

As a result, Eqs. (g]) and ([7]) become invalid in the vac- 
uum mass basis, and all the following derivations are in- 
valid. In practice, however, the adiabatic/precession so- 
lution may still be a good approximation even when 6 is 
large. This is because 9i3 ~ and the transformation 
matrix U is almost real. In fact, numerical simulations 
in Ref. ^2^ show that, at least for the parameters em- 
ployed in those simulations, varying the CP phase S has 
little effect on the final neutrino energy spectra except 
for changing the relative mixing of p and r neutrino fla- 
vors. The possible effects of CP violation in stellar col- 
lapse have been discussed in a different context (see, e.g., 
Ref. H). 

We also have demonstrated that the adia- 
batic/precession solution obtains even in the presence 
of a dominant matter background and a large mu-tau 
term. When the matter term is much larger than the 
vacuum mixing term, the presence of the ordinary 
matter only reshuffles the neutrino states in which the 
spectral swapping has the most dramatic manifestation. 
For the supernova environment, this means that the 
regime where collective neutrino oscillations occur is 
solely determined by the neutrino fluxes and does not 
depend sensitively on the matter density proflle. This 
is in agreement with the previous analysis in two-flavor 
scenarios jl^. The supernova neutrino signals observed 
on earth will depend of course on the matter proflle 
in the supernova. In part, this is because the MSW 
effect will modify the neutrino energy spectra subse- 
quent to the collective oscillations discussed here. This 
dependence of neutrino signals on the matter profile can 
provide important information on the conditions deep in 
the supernova envelope ^2^, iH, . 
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